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We discuss the scaling properties of free branched polymers. The scabng behaviour of the model is 
classified by the Hausdorff dimensions for the internal geometry : and dn, and for the external 
one : Dl and Dh- The dimensions dn and Dh characterize the behaviour for long distances while 
dh and Dl for short distances. We show that the internal Hausdorff dimension is cíl =2 for generic 
and scale-free trees, contrary to dn which is known be equal two for generic trees and to vary 
' between two and infinity for scale-free trees. We show that the external Hausdorff dimension Dh is 

directly related to the internal one as Dh = oidn, where a is the stability index of the embedding 
weights for the nearest-vertex interactions. The index is a = 2 for weights from the gaussian domain 
of attraction and O < a < 2 for those from the Lévy domain of attraction. If the dimension D of 
the target space is larger than Dh one finds Dl ~ Dh, or otherwise Dl = D. The latter result 
means that the fractal structure cannot develop in a target space which has too low dimension. 
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Introduction 



> 

\l ■ In recent years the theory of random geometry jl] has become a powerful method of investigating problems in 
' many areas of research ranging from the statistical theory of membranes |^ , branched polymersj |, ^ and complex 
networks Q to fundamental questions in string theory |l^ and quantum gravity |ll|, |l^, [l| ] . 

Those problems have in common that they can be described by a dynamically alternating geometry which undergoes 
"p-, fluctuations of a statistical or quantum nature. The dynamics of such fluctuations can be modelled using the concepts 
d of the statistical ensemble and the partition function in a similar way as it is done in particle physics by the methods 
^ ', of lattice field theory. 

I ' Contrary to lattice field theory where the partition functions run over field configurations on a rigid geometry, the 
■ O ] geometry itself is variable here. Since the geometry is dynamical many new features occur like for instance geometrical 
^ correlations or the influence of the random geometry on the fields living on it. 

^ ' Similarly as in field theory where the concepts of universality, critical exponents, correlations, etc are independent 
• • , of whether one discusses a field theoretical model of magnetism or a quantum theoretical model of particles, also 
. ^ in the theory of random geometry many questions are independent on details and may be addressed using generál 
k> methods. General concepts can be best developed on an analytically treatable model. In field theory the role of a 
j_] ■ test bed is played by the Ising model while in random geometry by the branched polymer model |l^, |l^, [l7| . 
', Despite its simplicity the branched-polymer model has a rich phase structure exhibiting different scaling properties 
of the fractal geometry and the correlation functions. 

The model has internal and external geometry sectors similarly as the Polyakov string ||]. As the Polyakov string, 
it can also be interpreted as a model for quantum objects embedded into a D-dimensional target space or a model 
of quantum gravity interacting with D-scalar fields. The term quantum gravity refers to an Euclidean Feynman 
integrál expressed by a sum over diagrams representing the nearest-neighbour relations between points of a discrete 
set. In more realistic models the sum runs over higher dimensional simplicial manifolds and can be interpreted as a 
regularized Feynman integrál over Riemannian structures on the manifold The intuition which one can 

gain from the analytic solution of the branched-polymer model is very helpful for considerations of more complicated 
models. In fact, the model has proven already many times to be extremly useful to test and develop various ideas 
concerning random geometry [pď|-|18 



In addition to this generál interest in this model there is a specific motivation which is related to the reduced super- 
symmetric Yang-Mills matrix model [p^ . This model was introduced as a non-perturbative definition for superstrings 
pO| and referred to afterwards as the IKKT model. The one-loop approximation of this model leads to a model of 
graphs which have as a back-bone a branched-polymer with power-law weights for the link lengths . The IKKT model 
is believed to provide a dynamical mechanism for the spontaneous breaking of the Lorentz symmetry from ten to 
four dimensions pl[ |. If one tries to understand the breaking in terms of the one-loop level approximation one finds 
it to be related to the fractal properties of the branched-polymers, which have the Hausdorff dimension equal to four 
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pl| , p2[ . The question of the spontanous symmetry breaking was investigated also by many other methods with help 
of which one was able to gain an insight into underlying physical mechanisms p5|-||3^. 

In the simplest case one considers generic trees with the nearest-vertex interactions given by gaussian weights 
or weights which belong to the gaussian universahty class jj, ||. In this páper we also discuss the nearest-vertex 
interactions, but we extend the discussion to non-gaussian weights in particular to power-law weights p^ . The 
energetical costs to generate long links on the polymer is then much smaller than for gaussian ones. In the extremal 
situation, when the power-law exponent a of the link length distribution ~ x^^^" lies in the interval a € (O, 2), very 
long links are spontanously generated on the tree and their presence shifts the model to a new universality class which 
can be called the class of Lévy branched polymers that similarly as Lévy paths exhibit a new scaling behaviour. 

The scaling properties and the universality class of the model depends also on the internal branching weights of the 
trees . Under a change of the weights the model may undergo a transition from the phase of elongated trees with 
the internal Hausdorff dimension dn ~ 2, known as generic trees, to the phase of collapsed trees with dn = oo which 
are localized around a singulár vertex of high-connectivity ^ . In betwcen there is a phase of scale-free trees 
which may have any Hausdorff dimension between dn = 2 and d// = oo |Q, |l^ . We show that the internal properties 
decouple from the embedding in the target space but on the othcr hand that they strongly affect the embedding : the 
external Hausdorff dimension is proportional to the internal one Dh = adu- 

In this páper we cover the entire classification of the universality classes of branched polymers with the nearest- 
vertex interactions. We hope it can be treated as a starting point in the discussion of models with more complicated 
interactions like for example those with correlations between ncighbouring links, or sclf-avoidance in the embedding 
space. 

Many pieces of this classification have been discussed for the gaussian branched-polymers already |l5|, ^ . 

Several well-known results have been summarized within the appendix in which we present a systematical treatmcnt 
of the internal geometry in terms of generating functions. 

The extension of this classification to weights with power-law tails and to the case when the internal geometry is 
non-generic, is presented in the main-stream of the text. Emphasis is put on calcuations of the two-point functions. 
Throughout the páper we also stress the factorization property of the internal and external geometry, which allows 
US to clearly separate the discussion of the internal geometry bcfore considering the entire model. It also permits us 
to reveal many interesting relations between the correlation functions of the external space to those for the internal 
geometry. 



The model 



We consider a canonical ensemble of trees embedded in a D-dimensional target space. The partition function of the 
ensemble is defined as a weighted sum over all labeled trees with N vertices. The set of such labeled trees, containing 
N^~^ elements, will be denoted by TJv- The statistical weight of a tree is given by a product of an internal weight 
Wt which depends only on the internal geometry of the tree, and an external one which depends on the positions 
of the (tree) vertices in the target space. We shall consider trees with nearest-neighbour interactions for which the 
partition function reads : 

1 r ^ 

ll^^^x., Hfir-,) . (1) 

■ TeTjv i=l (jk) 

The external weight of a tree is a product of link weights f{rjk) which depend exclusively on the link vector fjk —Xj~Xk- 
Saying alternatively, the energy cost of the embedding of the tree in the target space is a sum of energy costs of the 
independent embedding of links. The second product in (|^) runs over the set of (unoricnted) linked vertex pairs 
denoted by {jk). 

The most natural choice of the internal weights is Wt = 1. We could entirely stick to this choice of weights, but 
since we want to discuss the problém of universality we also want to check whether a modification of the weights will 



change the scaling properties and hence the universality |18 



Here we will restrict our considerations to internal weights which can be written as a product of weights Wq for the 
individual vertices : 

Wt = Wq^Wq^- ■ -Wq^ . (2) 



Each vertex weight only depends on the degree of the vertex that means the number of links emerging from it. The 
internal properties of the model are determined when the whole set of branching weights {wq} for w = 1, 2, . . . , cxd is 
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specified. We demand that : 

Wl>0, Wq>0, Ws>0 (3) 

for all q = 2, ... ,00 and at least one s > 2. lí wi were zero Wt would vanish for all tree graphs whilc if all Wq for 
q > 2 were zero then the weights Wt would vanish for all trees except chain-structures. 

Note that the model is invariant with respect to translations in the external space : Xi —i- Xi + 6. Because of 
the translational zero-mode the partition function (Q) is infinite. One can make it finite by dividing out the volume 
V ~ J d^x oí the translational zero-mode : 

ZN = ^. (4) 

This can for example be realized by fixing the position of the center of mass of the trees. 

Trees which can be obtained from each other by a permutation of the vertex labels contribute with the same 
statistical weight. For a tree with N vertices there are iV! such vertex permutations. In order to avoid overcountings 
one introduces the standard factor to the definition of the partition function (|l|). This factor divides out the 

volume of the permutation group of the vertex labels. The number of all labeled trees counted with this factor 
N^~-^/N\ ~ 7V~^/^e^ is exponentially bounded in the N ^ 00 limit. If one defines the grand-canonical partition 
function 

00 00 

N=2 N=2 

one can see that it is well defined as long as the fugacity g is smaller than the rádius of convergence of the series, which 
in the particular case Wt = 1 is equal go — e~'^" = e'^^. More generally, as long as zjv grows only exponentially for 
large N the grand-canonical partition function has a non-vanishing rádius of convergence and hence one can savely 
define Z. 

The statistical average of a quantity Q defined on the ensemble (|l|) is given by : 

11 r ^ 

(Q>^^— ^ E ll^^^^x, l[fif,k)Q. (6) 

^ • Ter„ 1=1 (jk) 

For translationally invariant quantities the averages are proportional to the volume V of the translational zero-mode. 
For such quantities one should rather speak of an average density per volume element of the target space : {Q)n/V 
which is a finite number. In particular (l)Ar/y = 1. 

We will frequently distinguish between the internal and external geometry of the trees. By the former we mean 
the connectivity of the corresponding tree graph, by the latter its embedding in the external space. For example, the 
internal (geodesic) distance between two vertices i and j of a graph is defined as the number of links of the shortest 
path connecting them, while the external distance is given by the length of the vector Xi — Xj. Note that the path 
between i and j is unique for tree graphs. Thus the length of this path, i.e. the number of its links, determines the 
internal geodesic distance. 

The properties of the embedding in the external space depend on the link weight function /(r) (see (^). We 
consider isotropic weights depending only on the link length. That means /(r) = /(r), where r = |ř^. We further 
assume that f(f) is a positive integrable function. Without loss of generality, we can then choose the normalization 
to be : J d^r f{r) = 1. This allows us to interprete /(r) as a probability density. 



Correlation functions 



The fundamental quantities which encode the Information about the statistical properties of the systém are the 
correlation functions. For the canonical ensemble with the partition function the m-point correlation functions 
are defined as : 

/ m AT \ 

\fc=l afc=l / jY 

where the brackets ( . )^ on the right hand side denote the average over the ensemble (^. If one multiplies all sums 
in the product in (R) one obtains a sum of terms being products of delta functions. The prime in eqviation (M) means 
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that all terms which contain two or more identical delta functions are skipped from this sum. This exclusion principle 
applies only to the situation when any two arguments of cj^^ (^Ai , • ■ ■ , ) are identical. 

What we will show now is that the problém of determining the m-point correlation function can be devided into two 
sub-problems. The first step is to determine the internal two-point function. This can in generál be done independently 
of a particular choice of the embedding weights f(f). The second is to use the Information encodcd in the internal 
two-point function to determine the external properties of the trees. 

In order to see that the internal geometry of the trees does not depend on the choice of the link weight function, 
consider a tree graph and calculate the foUowing two integrals for this tree : 

N 



í n d''^^ n / (^j'^) • '^(^- - ^-4) = 1 , (8) 

^=l (jk) 

N 

lld^^x, llfir^k) ■ S{xa-XA) S{xb-XB) - UXb - Xa). (9) 

í=i ijk) 

The first integrál corresponds to the embedding weight factor for a tree whose a-th vertex is fixed at the position Xa 
in the target space. Since the model is translationally invariant the result of the integration does not depend on the 
position. This result can be obtained by changing the integration variables from position vectors of all vertices i ^ a 
to link vectors rjk — Xj — Xk for which the integration completely factorizes. The Jacobian for such a change of the 
integration variables is equal one. 

The second integrál gives the weight factor for a tree whose vertices a and b are fixed at the positions Xa and 
Xb in the external space. The result depends only on the difference X = Xb — Xa and the number of links, n, of 
the path connecting the vertices a and b. If one now changes the integration variables from vertex positions to link 
vectors, as before, one can see that all integrations, except those for links on the path between a and 6, factorize. The 
sum of the link vectors on the path is restricted to X = Xb — Xa- If we label the links of the path by consecutive 
numbers from 1 to n, we can write : 



fn{X) 



/n / \ 

X{{á^nf{n))5[Y,ra-x\ 
1=1 \a=l / 



d^D 1" -i? 

^ í{p) e^^f^ , (10) 



(27r)^ . . 

where /(p) is the characteristic function of the probability distribution /(r) : 

m f d^r f{r) = {é^f . (11) 



It is important that the results of the integrations (^) and (^ do not depend on the internal geometry of the undcrlying 
tree graph. In particular, using (||) and (||) we find the partition function (Q), the one-point and two-point correlation 
functions to reduce to the foUowing form : 

= W^T.WT. (12) 
G^^X^) = 1, (13) 

N-l 

G^^\Xa,Xb) = G^^\X) = /"(^) 9N\n) • (14) 

n=l 

(2) 

We have denoted the (canonical) two-point function of the internal geometry by g)^ (n) : 

9'N\n) = E ]^ E '^la-.Kn I . (15) 

■ TeTjv \ Q,ber 

It is normalized to unity : X]í^=o^ ^ 1- ^^'^ normalizcd internal two-point function gives us the probability that 

two randomly chosen vertices on a random tree of size TV are separated by n links. In the last formula the geodesic 
internal distance between the vertices a and b is denoted by |a — 6|. 
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The expressions ( p^ for zjv, ( p^ for G^^\x) and ( p^ ) for g^^\n) are independent of the link weight factor /(r). 
Thus the propcrties of the internal geometry, as mcntioned, can be considered independently of and prior to the 
embedding. 

The external two-point function Gj^ {X) can be interpreted as the probabihty density for two random vertices on 
a tree of size N to be embedded in the external space with the relative position X = Xb ^ ^a- The probabihty 
normahzation condition reads : 

" d^XBG^^\XA,XB) = G^^\XA) = l. (16) 

The right hand side of (|lj) can be understood as a conditional probabihty. First, we choose two random vertices on 

a tree of size N and calculate the probabihty g)^ {n) that there are n hnks on the path connecting thcm. For this 
path, which is a random path in the embedding space consisting of n hnks, we can calculate the probability (density) 
fn{X) that its endpoints are located with the relative position X ~ Xb — Xa- Since the internal geometry decouples 
from the external one, the probability densities fn{X) and g)^ (n) are independent of each other and can hence be 
calculated separately. 

Similarly, higher correlation functions can be obtained from the corresponding internal correlation functions. For 
example, the three-point correlation function is : 

G%\Xa,Xb.Xc) = 9NÍ^a,m,n,)x (17) 

ria jí^ifa ,í^c 

X J d^X [Xa - X)U, {Xb - X)U^ (Xc - X) . 

The three paths ab, bc and ac between the vertices a, 6 and c of the tree can be decomposed into three pieces, 
námely am, 6m, cm between them and the common middle point m. The summation indices Ua, rib, ric denote the 
internal lengths of these pieces, and X denotes the position of the common vcrtcx in the external space. The internal 
three-point function then reads : 



gfina.m.n,) = —-i- V PFt x (18) 



/ j '^la— 6| ,na+íib ^1 ř?— c| ,?ib+íic ^|a^c| .na+Tic 
a,fc,ceT 

One could extend this construction further. 

Note that the most important piece of Information is already encoded in the two-point function and is inherited by 
the higher correlation functions. In fact, one can directly derive the higher correlation functions from the two-point 
function, using a simple composition rule for the tree graphs which enormously simplifies in the grand canonical 
ensemble. In the next section we shall thus concentrate on the two-point function. 



Fractal geometry 

f2) Í2) ^ 

The canonical two-point correlation functions g]^ (n) and G ]^ {X) contain the Information about the fractal struc- 
ture of the internal and external geometry, respectively. The average distance for the internal geometry, given by the 
average number of links between two vertices on the tree, is the first moment of the probability distribution (na) : 

n 

One expects the foUowing scaling behaviour for large TV : 

{n)N-X^I'^«. (20) 

The exponent dn relates the systems average (internal) extent (n) to its size X and is thus called the internal 
Hausdorff dimension. This exponent controls the behaviour for large distances growing with the systém size X . One 
can also introduce a local definition of the fractal dimension for distances in the scaling window 1 ^ n ^ jVi/^^ff. 
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The scaling window contains distances between the scale of the ultraviolet cut-off and below the infrared scale set by 
the systém size. This is a sort of thermodynamic definition which becomes valid locally for sufficently large N. In a 
large systém one can be interested in how the volume of a local ball (or sphere) depends on its rádius n. The volume 
of the sphere can be calculated as the number of vertices lying n links apart from a given vertex : 

g^l^^(n) - rA-i for 1< n < A^^/''" . (21) 

The definition d^, is more practical for a local observer, for example, someone who lives in a fractal geometry and 
wants to determine its fractal dimension. The globál definition dn is accessible only for an observer who can survey 
the whole systém from outside Q . 

In a similar way we can dchne the external Hausdorff dimension. In order to do this we first have to introduce a 
measure of the systems extent in the external space. Such a measure is provided by the gyration rádius : 

^' = ^J2i^^- = ^ E(^» - ^CMÝ , (22) 

where Xcm — TIí-^í/-^ target space position of the systems center of mass. The statistical average of the 

gyration rádius is directly related to the two-point function, námely : 



^{R')n^ J d^X X' G'j^>{X). (23) 

Since the gyration rádius is a translationally invariant quantity we have to normalize it with the total volume of the 
target space and rather refer to its average density. The external Hausdorff dimension Dh can then be read off from 
the large N behaviour : 

y/WU - ■ (24) 

The dependence on the volume V is hidden in an A^-independent constant which is not displayed in the last formula. 
The symbol ~ referes to the leading behaviour. 

For trees of sufBcently large size N one can also definc a local fractal dimension Dl of the external geometry 
by measuring the average number of vertices within a spherical shell of rádius X = |X| from the scaling window 
Xuv <^ X <ti aN^^-^" , which is defined above the ultraviolet cut-off scale and below the infrared scale. As foUows 
from the definition (|^) for the case of two-point function, the number of vertices within a spherical shell of width dX 
is given by the two-point function [ p2[ : 

n{X) dX - X^-^ {X) dX - X^^-i dX . (25) 

Here we have used the fact that the the two-point function is spherically symmetric, i.e. G)^ {X) = G)^f {X). 
The integrál over the I?-dimensional angular part is included in the proportionality constant. In the large N limit 
one expects the existence of a window Xuv ^ X <^ aN^^^" where the two-point function exhibits the scaling 
behaviour |Q : 

G^^\X)^X-'. (26) 
If S is negative then G)^ {X) behaves as a slow varying function of X, which can in a narrow range and with some 

(2) 

corrections be viewed as a constant G)^ (X) ^ 1. Thus depending on the value of 6 we have : 

In the first case (S < 0) the number of points in a spherical shell n{X) grows with the power of the canonical dimension 
D. Only in the second one {6 > 0) the fractal nature leaves traces in the calculation of D^. 



Universality classes and singularity types 

In this section we will shortly summarize results concerning the classification of the scaling behaviour according to 
the internal geometry of the trees [0, |l8) . One defines the critical exponent 7 of the grand-canonical susceptibility 
via : 

= - A^-- , (28) 
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where A/x = fj, — fj,Q controls the behaviour of the partition function at the rádius of convergence Qq = e~^° of the 
series (||). Here /íq is the critical value of the chemical potential. If 7 is positive the susceptibility Xfj. itself is divergent 
at /íq. If it is negative, the right-hand-side of the last equation should be understood as the most singulár part of the 
susceptibiUty, which after taking higher derivatives, will give the leading divergence. The primary classification of the 
universahty classes for models of branched geometry is based on the value of the susceptibility exponent 7. 

The susceptibility exponent gives the subexponential behaviour of the canonical coefRcents zn for large N : zm ~ 
iV^~^ exp(/io-/V)- Indeed, if one inserts this form to the definition of the partition function (|5|), one obtains : 

= S = E N^^.NC-^'' - E iV^-'e-^^^ ^ A^-^ . (29) 

^ N N 

The susceptibility is proportional to the hrst derivative of the grand-canonical partition function $ for planted 
rooted trees, defined by equation (|77| ) in the appendix : ~ d^/du. The reason, why this relation is useful is, that 
there exists a closed relation - a so-called self-consistency relation (pOh - for $ : 

which can be inverted for í> = $(5) and from which one can extract the singulár part of $ : ^ ^ ^g^~^ and hence 
also of X/i- Note that the dcnominator of ( ^0| ) is nothing but the first derivative of the potential V[^) defined ( |7^ ) 
within the appendix. 

One can invert the function .g($) in the region where it is strictly monotonous. Generically .g($) grows monotonically 
from zero for $ = O to some critical value at $0 where 5'($o) = 0. Clearly the inverse function <i> has a square 
root singularity at 50 : $ then. It follows that 7 = 1/2 for the class of generic trees. 

The region of the monotonous growths of the function on the right hand side of (|3^) may be limited by $0 being 
the rádius of convergence of the series in the dcnominator of 5(*í'). 

The inverse function $(g) is then singulár at — g{^o) with a singularity inherited from the singulár behaviour 
of g{^) at $0- It can be shown that in this case the susceptibility exponent is negative and the corresponding trees 
are collapsed. 

In the marginal situation the two conditions which limit the region of the monotonous growth of g{^) work 
collectively at a point $0 being at the same time the rádius of convergence of the series in the dcnominator of (?(<ř) 
and the zero of the derivative g'{^o) = 0. In this case the exponent 7 can assume any value within the interval 
[0, 1/2). Trees which belong to this class are called scale-free |7|. 

The three classes correspond to different scaling behaviours of the two-point function as will be discussed in the 
foUowing section. 

Internal two-point function 

In this section we shall calculate the two-point function for the internal geometry. This function will enable us to 
determine the scaling and the fractal properties of the internal geometry of tree graphs. As we will show they are 
different for generic, collapsed and scale-free trees. 

In the appendix we deduce an explicit formula ( p7| ) for the grand-canonical two-point function í/*^^^ (/i, n). This 
function is singulár for A/i = /i — /íq 0+, and singularity is related to the large iV-behaviour of the canonical two- 
point function ÍJ^-* (n). The singularity of G^'^^{p,n) can be determined directly from the identity ( p7| ) by inserting 
the most singulár part of A<í> to V'{^) and V"{^). Here we will show an alternativě way using a standard scaling 
argument from statistical mechanics pq . Denote the singularity exponent of the two-point funtion by v : 

^(2) (^, n) ^ exp (-c(n + í)An'') , (31) 

where c is a constant which only depends on the particular choice of the weights Wq. The exponent v is usually called 
mass exponent. Summing over distances n we obtain the susceptibility ( ^ ) : 

Xa. = E ^^'^ í'^' ") ~ / ^" (-cnA/i"-) ~ A^^-'' . (32) 

According to the definition (p9|), the susceptibility exponent is 7. Thus we have 



= 7 . 



(33) 
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This relation is the Fisher scaling relation for this case. Since we have already determined 7 we do not have to calculate 
v additionally. The scaling argument given above holds only for positive 7, since in this case the susceptibihty is 
divergent and the divergent part dominates the small A/z behaviour. For negative 7 the left hand side of (|2|) is not 
divergent which means that it behaves as a constant as A/i goes to zero. In this case, if one compares the resuh of 
the integration ( ^2| ) to the leading behaviour of (p^, one shaU effective see that v — Q. This is what happens in the 
collapsed phase. 



Now, inserting the most singulár part of the grand-canonical two-point function (31 ) to the inverse Laplace transform 
(p^) we can deduce the large N behaviour of the canonical two-point function : 



where 



g^^\n)^e+^'>^L,{cn,N), (34) 

LAcn,N)^— dCe-™í"+«^ (35) 

27^^ J^^_,oo 



is the Lévy distribution with the index v, the maximal asymmetry and the range C = cn cos (71W2) 
The large N asymptotic behaviour of L^^cn, N) with v < 1 is given by the foUowing series ||24| : 



^ ' ^ ttN ' vTv^y r(i + fc) ^ ^ ' 



k=l 



For large N and fixed n, the first term dominates the behaviour of the series : 

, . w^oo lyTjiy) sin(7ri/) cn 

LAcn,N) ^ (37) 

We see from the formulae ( ^^ and (|35| ) that the two-point correlation function in the large N limit is effectively a 
function of the argument cn/N^ . Indeed, if one changes the integration variable ^ in (^) to ^' = ^A'' one obtains : 
L^{cn,N) = N~^L^{cn/N'' , 1) = N^^l^{u), where l,y{u) is a function of a single argument u = vcn/N'^ . For later 
convenience we also included v into the definition of the universal argument. Using the saddle point approximation 



to the integrál (35) one can find that for large u = vcnjN^ the function ljj(u) leads to: 



H - j^e+^^^^Uu) = ^"'^ exp {-hu'^) , (38) 

where a~ 1/(1 — z^) and 6 = (1 — v)lv- The average internal distance between two vertices can then be calculated : 

{n)N = 2_^ ngy (n) = .00 , ... ■ (39) 

„ Í^C du l^u) 

Comparing the Aí^-dependence on the right hand side of this equation to the iV-dependence on the right hand side of 
equation (^0|), which defines the internal Hausdorff dimension dn, we see that dn is the inverse of v : 

dH = - = -. (40) 
u 7 

Thus, the Hausdorff dimension is dn — '2 for generic trees. For scale-free trees dn changes continuously from 2 to 
00 smce 7 belongs to the interval [OJ/2) then [5[ |l|. 

On the other hand we see from ( |37| ) that for large TV and small n the normalized two-point function grows linearly 
with n, i.e. : 

9^Nl^{n)-n. (41) 

The normalization coefhcent behaves as c/N in the large N limit. Since the sum over this function is proportional 
to the number of vertices in the distance n from a given vertex, the last formula tells us that the local Hausdorff 
dimension is cŽl = 2. We see that locally for sufficently large N, it is difficult to distinguish the scale-free trees from 
the generic ones by measuring short internal distances, since both classes have the same internal Hausdorff dimension 
d^ — 2. One has to go to large distances to see different scaling properties depending on the type of the ensemble 
(pq). For collapsed trees the Hausdorff dimension is infinite. In this case, the universal scaling argument u of the 
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two-point function is proportional to n but it does not depend on . This is related to the fact discussed before that 
the efFective value of the exponent v is equal zero. 

The saddle point approximation (38) actually gives the exact result for ly — 1/2 for the whole range of u. The 
reason for this is that in this case the integrand of the approximated expression (|35| ) is gaussian. For some specific 
values of ly one can express the Laplace transform (pq) in terms of speciál functions. For example for v = 1/3 : 



(42) 



where u — vcn/N'^ . For large u the saddle point formula (38) coincides with this one, while for small u the two 
functions deviate a littlc from cach other. 



Gaussian trees 

Now we can determine the properties of the external geometry of gaussian trees. In this case, the weight f{x) in 
the partition function (|^) for embedding a link x is given by a gaussian function. The function has a vanishing mean : 

/(f) = (27ra2)-^/2exp(^~^) . (43) 

In other words, for gaussian trees the link vectors are independent identically distributed gaussian random variables. 

As a consequence, the probability density fn{X) ( p^ ) for the endpoints of the path of length n on a tree to have 
the relative coordinate X is given by : 

. (2™.',-'^ exp . (-L)% . ,44, 

This foUows from the stability of the gaussian distribution with respect to the convolution. Inserting the function 
/„(X) to the formulae (H), ( pŤ|), etc. we can determine the multi-point correlation functions for gaussian trees. In 
particular, if we insert ( ]4^ to (|14|), we obtain in the large N limit : 

Here we used the same approximation for the internal two-point function as in the discussion of equation ( p4[ ) in 
the appendix. This is a good approximation for large N . Additionally, we substituted the upper limit — 1 of the 
summation over n by infinity. This introduces small corrections which disappear exponentially in the large N limit. 

In order to measure the external Hausdorff dimension Dh we have to determine the dependence of the expectation 
value {R'^)n of the gyration rádius on the systém size N. The expectation value can be calculated by integrating the 
two-point function over X as in equation (p3|). If one first integrates over X before summing over n one obtains : 

'^iR')N-'-^y:n'e^p{---] . (46) 



V ' N 



One can approximate the right-hand side by replacing the summation from 1 to cx3 through an integration over the 
whole positive real axis : 

rAR')N-^ dyy^e^p ~JL\^Da\ -. (47) 



V' N Jn " " V 2N J V c 



We see that the typical extent of the systém \/{ĚF)n grows as iV^/* and hence the Hausdorff dimension for generic 
gaussian trees is Dh = 4. 

More generally, in order to determine the dependence of the expectation value of the gyration rádius on the systém 
size for any type of trees one can first calculate the second moment of the function /„(X) : 

(X2)„ = / d^X X^UX) , (48) 
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which corresponds to the average extent of the path built out of n links of the tree. The insertion of this result to 
(H) yields : 

l(i?2)^^^(X2)„g(^2)(n). (49) 

n 

Since for gaussian weights the second moment (X^)^ is proportional to n, i.e. {X^)n ^ the foUowing relation holds 

M ■■ 

{R')n ^ ^ng(,''(n) ^ (n) ^ ^ N^'''- , (50) 

n 

from which we conclude that the external and internal HausdorfF dimensions are related by : 

Dh = 2dH (51) 

for gaussian trees. This relation holds for generic, scale-free and coUapsed trees. This for example means that the 
Hausdorff dimension Dh of coUapsed trees is inflnite, or in other words, that the target space extent of the systém 
does not change with the number of vertices on the tree. 

Let US come back to generic gaussian trees. We will calculate the local Hausdorff dimension Dj^ and compare it with 
Dh — 4. The starting point of this calculation is equation ( p7| ) which relates the number of points within the shell of 
rádius between X and X + áX to the behaviour of the two-point function in the scaling window Xjjv <C ^ <C ^/_r, 
where Xjjy cr, Xjr ^ aN^^^. 

The two-point function ( |45| ) is a decreasing function. It has a cut-off sá X ^ aN^/^ as follows from the scaling 
arguments. The large n-part of the sum (^) over n : n 3> 1 can be approximated by an integrál over n. This part of 
the sum has a significant contribution ifC = X^/(T^^l. Thus for X ^ cr the dominant dependence of the sum on 
X can be approximated by : 

l'C2\ÍŇ 

G^^ {X) ^ án ni-^/^e-í/" ^ C'^""/" ■ (52) 

J Cl 

with some constants ci, C2. The uppcr limit of the integrál comes from the term exp(— C7i^/2A^). The cxact shapc of the 
integrand at large n ^ \fŇ is unimportant for D > 4, because the dominating behaviour G^^\x) ^ ^2-d/2 ^ X"^^^ 

(2) 

is due to the lower limit of the integration. This scaling form of G)y {X) breaks down for short distances X of order 
a and for large X of the order of the infrared cut-off aN^^'^. When X is of the order a the integrand is a sum of 

(2) 

gaussians of widths larger than X and hence G)^ {X) is a slow varying function. 

For D < 4 the regime changes. The divergence at small n disappears and the tcrms for large n, n ^ \ÍŇ, dominate 
in the sum. The sum (^), viewed as a function of X , looks as a sum of gaussians whose arguments X are maximally 
of the order of the widths. This is a slow varying function of X for X^/ct^a/ŤV- Hence we expect that as long as 
X ^ N^/^ it is almost constant g''^\x) ^ 1. As an example we performed the sum ( ^ ) numerically for N up to 

10^. The results presented in the figures |^ and || corroborate the anticipated behaviour of Gj^"* {X) by the arguments 
given above p0| . 

As a consequence we see that the number of points within the spherical shell (^^ depends on the rádius X as : 

r X^-i for D<4 
^{x^ íor D>A. (5^) 

This leads to the foUowing result for the Hausdorff dimension : 

^ f D for !?< 4, 

^^ = (4 for D>4: (54) 

In other words, the fractal dimension Dl measured by the local observer, is equal to the globál one, i.e. Dh — Dl, if 
the dimension of the target space is large enough. If the target space dimensionality is too small, the fractal structure 
cannot develop. One can understand this in the foUowing way. For trees embedded in a D < 4 dimensional target 
space, vertices of the tree lie in a ball with a rádius proportional to iV^/"*. There are N vertices within the ball while 
the volume of the ball is proportional to 7V^/^. This means that for large N, the vertices deep inside the ball are 
densily and uniformly packed. A local observer who surveys only a small region far from the ball boundary will see 
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FIG. 1: The normalized two-point function for D = 2 and for A'^ = 10^, 10*, 10^ and 10®, from left to right, respectively. The 
functions are constant in the region of small X. This region extends to some cut-off whose position grows with the power N^^*. 




X/a 



FIG. 2: The normahzcd two-point function for D = 10 and for N — 10"^, lO**, 10'' and 10^, from left to right, respectively. The 
functions are constant in the region of small X. Thcn thcy dcvolop a scaling part (lincar in the figuře) in which they behave 
as cx X~^. For a comparison of the slopes we also displayed (solid line) a pure power-law ex X~^. 



uniformly distributed vertices in a D dimensional spacc. As a conscqucncc hc or shc will mcasuro Dl — D. The 
situation changes for _D > 4, because then the volume of the ball is proportional to N^/^ and hence grows much 
faster than the the number of the vertices. In the large N limit the volume of the ball will therefore be large enough 
to let the systém develop a loose fractal structure. 

Similarly, we expect that for the scale-free trees, the local Hausdorff dimension is 

j D ÍOT D<Dh, 
^^ = \Dh for D>Dh. 

Since the Hausdorff dimension Dh is infinite for collapsed trees, one cannot deíine a local Hausdorff dimension Dl 
in the same manner as above, because the infrared cut-off is a constant in this case. 
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Lévy trees 



So far we have only considered gaussian embedding weights f{x) (p^). In this case links typically have the length 
a and one can hardly find a link on a tree longer than Su. In other words the energy costs for the embedding of long 
Unks are so high that those links do not appear. One can, however, consider models with weights which allow long 
links. In the next section we will discuss this issue in a more generál context, while in this section, as toy models, we 
will consider models of trees embedded in D = 1 dimensional target space with the weigths given by a symmetric Lévy 
distribution |2^ . Despite their simplicity, the models with those weights already basically capture all interesting 
features of more complicated models. The weights read : 

/+00 
dp exp(-A"|pr-ipx) , (56) 
-oo 

with a from thc interval (O, 2]. In thc limiting case a = 2 is a gaussian distribution with thc width a = \/2A, and for 
a = 1 is a Cauchy distribution. 

The weights are symmetric stablc distributions with the stability index a. Herc we are intcrestcd only in symmetric 
functions f{x) — /(—a;) because the links are unoriented. This implies that any function defined on them has the 
property : f{r,j) = lLr'^ij) = /(-"^jO = fi^jz)- 

The distribution (^6|) is stable with respect to the convolution : 

r+oo 

La{A,x)^ dxi La{Ai,Xi)La{A2,X - Xi) , (57) 



where A" = Af + A2. If one repeats this for the convolution of n identical terms to calculate fn{X) ( |1 1| ) one can see 
that the function fn{X) is given by a rescaled version of the function for a single link f{X) : 

1 „ í X 



Now we can combine this scaling with the scaling of the internal two-point function which, as we know, is a function 
N~^l^{v) of a scaling variable v = n/N^ , to deduce the scaling of thc cxtcrnal two-point function ( p^ : 

cfiX) ^Y.í'.mň\n) --U^f (^) (^) . (59) 

The result of thc integration can be written as a function of an argument X/N''^" with some prefactor depending on 
N. As a consequence one expects the external Hausdorff dimension to be : 

Dh = a/v ^ adn ■ (60) 

The case a = 2 was discussed before. Despite similarities, the case a < 2 is different from the gaussian one, since in 
this case the distribution f{x) has a fat tail for large x : 

dx fix) ^ , (61) 

X x°' 

which according to the scaling ( |58| ) is equally important in fn{X) for any n . The second moment {X'^)n of the 
distribution f„{X) is infinite. As a consequence, also the gyration rádius is infinitc. One has to find an alternativě 
measure of the linear systém extcnt in order to define the Hausdorff dimension D^- A natural candidate for such a 
quantity is : 

^' = ]^El^'-^^l'' (62) 

id 

for q < a. The Hausdorff dimension can now be calculated from the large N behaviour of this quantity : 

{R'^)n - N^í'^" ■ (63) 
Using the same arguments as for the gaussian case, one can check that the foUowing relations hold 

^{R'^)n ^ jáX G^^(X) = 5](|Xr)„ gf{n) , (64) 
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where 

(|X|«)„ ^ JdX W - . (65) 

It foUows, that (i?9)Ar - TV-í/adff and hence Dh = adn, as already mentioned. 

Other trees 

We will continue the discussion of the one-dimensional case, i.e. D = í. The extension to higher dimensions D 
shaU afterwards be straightforward. The embedding weights for hnks may in generál be given by any normahzable 
non-negative symmetric function : /(x) = f{—x) > O such that : J dx f{x) — 1. 

We are interested in the emergence of the scahng properties for large N . From the considerations of the internal 
geometry we know that the internal distance between two random vertices on the tree, n ^ N^^^" , grows with A'^ 
unless the trees are collapsed. 

We also know that between those two random vertices we can draw a unique path on the tree. This path can be 
treated as a random path of n links. So, in a sense, we are interested in the probability distribution that the remote 
ends of the random path with n links have the the relative position X in the embedding space. In particular we 
are interested in the limit n — > oo. This probability distribution is given by fn{X). For large n the function fn{X) 
can be determined from the centrál limit theorem. Roughly speaking, if the second moment of f{x) exists, fn{X) 
approaches a gaussian distribution with the variance (t„ — y^a, otherwise fn{X) approaches the Lévy distribution 
( pq ) with the scale parameter A„ = n^/°'A. Thus, if a distribution has a power-tail f{x) ~ x~^~°' for large x, the 
limiting distribution fn{X) for large n will approach the gaussian distribution if a > 2, or the Lévy distribution ( [S^ ) 
if a < 2 The limiting case f{x) ~ x~^ belongs to the gaussian domain of attraction but it has a logarithmic 

anomaly of the variance which in this case does not grow as ^/n but faster, i.e. with some additional logarithmical 
factor of n. 

For a > 2 the approach of ,fn{X) to the gaussian distribution for large n is non-uniform and takés pláce in the 
centrál region of the distribution fn{X) for \X\ < X^, where X^ scales with n as 



X^ ^ h\/n logn . (66) 

Here h is some constant representing a scale of the distribution. Outside this region /n(X) deviates severe from 
the normál law, and in particular it preserves the power-law tail for X 3> ^* : 

dX/„(X)^^!^, (67) 

with a tail amplitudě proportional to n. In other words, for any finite n the power-law tail is present in the distribution 
fn{X). Therefore all absolute moments of order Q > a of this distribution are infinite : 

{\X\Q)„^oo (68) 

for finite n, and as a consequence also 

{RQ)n= j áX xQG%\x) = Y,{\X\Q)J^\n)^<^. (69) 

For N oo the sum is dominated by terms of large n ~ JSf^^-^" . For n ^ oo the distribution fn{X) becomes normál 
in the whole region from — oo to +00. Indeed, the ends of the centrál region X^, move to infinity faster than the 
variance cr„ ~ ^Jn and the contribution coming from the outside of the centrál region \X\ > \X^,\ disappears as : 

áX UX) \ . (70) 

\x\>x, n°'i'^ log ' n 

Thus the non-gaussian part including the tails becomes marginal and can be neglected. What is left over for n = 00 
is a gaussian distribution with all moments defined. For example, even integer moments Q = 2K are : (X^^)„ = 
{2K — 1)!! [na^]^ . Thus, after taking this limit the trees behave like gaussian ones. This limit is subtle, because as 
long as N (and n) is large but finite higher moments, Q > a, of the distribution are infinite. 
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Now we shall shortly discuss the model in D > 1 dimensions. As before we consider spherically symmetric distri- 
butions f{x) — f{x), X — \x\ which have a power-law dependence for large lengths of the link vectors f{x) ^ a;^^^" : 

á^x f{x) = d^D dx x°-^f{x) - dnodx x-^-" , (71) 

where dfio is the angular part of the measure. The main difFerence to the one dimensional case is that the effective 
power of the hnk length distribution changes by -D — 1 due to the angular measure x^~^ . Otherwise, the dependence 
of the scaling on a goes in parallel to the one dimensional case that is the distribution belongs to the gaussian domain 
of attraction if a > 2 and to the Lévy one if a < 2. The characteristic function ( pŤ| ) of the corresponding limiting 
distribution is spherically symmetric : exp(— where p = |p| and a = min{2,a}. For example, a distribution 
which has a power-law tail f{x) ~ x~^~^ belongs to the domain of attraction of the Cauchy distribution : 

f(^) = 7^ /d^? ^ im (72) 

A new effcct which arises in higher dimensions is a possibility of a spontaneous breaking of the rotational symmetry. 

(2) ^ 

The limiting distributions for N ~t oo and the two-point function Gjy {X) are spherically symmetric but the config- 
urations which contribute to them are not. The effcct is strong when a < 1 and is known from the considerations of 
Lévy paths |2^. If we have such an ensemble of N links, one can find a link whose length is iV^/" larger than the 
sum of the lengths of the remaining links. This link makes the systém look like a one dimensional systém, since the 
extent of the systém in the direction of this link is significantly larger than in the other directions. The effect becomes 
weaker when a is larger than one. Actually it is then seen for configurations which come from the large X-tail of the 

(2) 

two-point function G}^ {X). This configurations become marginal for a > 2 in the large N limit. However, as long 

(2) 

as N is finite, the probability of large X in Gj^ {X) is finitc and it strongly influcnccs the mcasurcmcnts of higher 
moments {R'^)n of the systém extent. In other words the higher Q is, the stronger is the contribution from the large 
X part of the two-point function and the more contribute systems, which are elongated, to this quantity. In the 
limiting case . Q a the main contribution to the moments {R'^)n comes from one-dimensional configurations. 

As mentioned, the branched polymer model with power-law weight arises as the one-loop approximation of the 
reduced super-symmetric Yang-Mills model pl[ . In particular for D = A dimensions, the embedding weights f{x) 
behave as f{x) ~ x~^~°' ~ x~^ for large link lengths x. This is the marginal case a = 2 which belongs to the 
gaussian domain of attraction. This means in particular that if one first takés the limit N oo then the gaussian 
branched polymers emerge for which all correlation functions are well dcfined. On the other hand if one determines 
higher correlation functions before one takés N oo one shall see they are divergent. 

In numerical simulations of the full matrix model one also observes power-law tails in the distribution of the systém 
extent and one dimensional configurations ||2^]-|2^. It is possible but not yet answered, that there also exists a 
Gaussian limit at large A'^ in this model. 



Discussion 



We investigated the model of trees embedded frčely in a D dimensional target space. We classified the scaling 
properties of the model by determining the fractal dimensions for internal and external geometry for the ensembles 
of generic and exotic trees including those which have fat tails in the distributions of branching orders and of link 
lenghts in the embedding space. We showed that, for freely embedded trees, internal geometry is indepedent of the 
embedding as a result of the factorization (|l^). On the other hand external geometry strongly depends on internal one : 
in particular the Hausdorff dimension for external geometry is proportional to that for internal one Dh = adn- The 
proportionality coefficent is given by the stability index of the embedding weights. For gaussian trees, in particular, 
it is equal two. We pointed out that the finite effects related to the presence of fat-tails lead to singularities of higher 
order correlation functions before the inifinite large A-limit is taken. This is a similar effect to the one observed in 
the IKKT matrix model. 

The branched polymer model captures many features of the more complicated models of random geometry. Despite 
its simplicity the model has a rich phasc structure : a generic phase of gaussian trees which have the Hausdorff 
dimensions Dh — 4, the phase of short trees with Dh > 4 coming from the embedding of scale-free and crumpled 
tree graphs, and the the phase with elongated Lévy branches, which has the Hausdorff dimension Dh < 4. 

Due to the simplicity, and the full control of the free case, the model of branched polymers which we discussed here, 
provides a good starting point for modelling effects of non-trivial embedding like those related to excluded volume 
effects or external curvature. Such effects violate the factorization introducing correlations between the internal and 
external geometry. A sort of back-coupling occurs. The external geometry affects the internal one, which modifies 
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and influences back the external one. For example, self-avoidance disfavours crumpled trees and hence changes the 
internal Hausdorff dimension which in turn will change the external one. 
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Appendix 

This appendix summarizes important relations of the internal geometry of branched polymer models ^ |8|, ^ . 
It is intended to make this article more self-contained. As we already mentioned, the part of the model related to 
the internal geometry decouples from the problém of the embedding and can hence be solved independently. It is 
convenient to introduce several generating functions to ease the calculations. Although many of the considerations 
made here are well-known, we deduce several relations for the generating functions which play important roles in the 
description of branched polymer models. Graphical representations of the generating functions turn out to be effective 
tools for the just mentioned deductions. Finally we will be able to calculate the partition function z^r ([l^ ) and the 

two-point function g^^\n) (|l5|) of the internal geometry. Note that we already discussed the scaling properties and 
the universality resulting from those calculations in the previous sections. 

Kecali that the internal properties of a branched polymer model (Q) only depend on the internal weight function 
Wt- For internal weights of the form (0) these properties are entirely determined when the wholc set of branching 
weights {wq}, obeying the conditions is given. The Information about the entire set {wq} of branching weights 



can alternatively be encoded in a single function of one real variable, námely (73) 



9=1 ^' 

As we shall show below, the scaling properties of the internal geometry are directly related to the analytic properties 
of this function. We will often refer to 1/ as a potential, since the most important generating functions can be written 
as derivates of V . 

In the first section we defined the generating function for the canonical partition functions to be : 

oo oo 



^ ZMg"" = 5J z^e-^^ , (74) 



which is nothing else but the grand-canonical partition function for the ensemble of trees with unrestricted size. One 
reason why it is convenient to introduce the generating function Z is that one can write a closed self-consistency 
equation for a first derivative of it, as we shall discuss below. The grand canonical partition function can be written 
as : 

oo 

^ = E M ^ ^^^^^"^^ ■ ■ ■ (^"'iv-i)^"-^ , (75) 

where Nq denotes the number of vertices of order q and the sum begins with N = 2 being the number of vertices 
on the smallest tree. Note that each vertex introduces a factor (gwq) to the (internal) weight of the tree in the 
grand-canonical ensemble. There are two derivatives of the generating function Z which will be useful : 

o ^ oo 
^ ^ N=l 

where 

dzN+i 1 NiWt 



1 ^ N,Wt 



dwi (iV + l)!^|^^^ wi 
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Clearly, Z^^^^ is a generating function for the canonical partition functions z^j^ = Nzn of trees with TV vertices that 
have one marked vertex. Intuitively the factor N in the sum can be viewed as a factor which counts the possible 
choices of marking one vertex on a tree with N vertices. The derivative $ is a generating function for the partition 
functions ipN oí branchcd polymers of size N having one (not counted) additional marked vertex of order one. We 
wiU refer to those uncounted vertices as external Hnes. Because we do not count the empty Unes in (77) the sum 
starts with = 1. 

We now introduce a graphical notation for the generating functions with the foUowing conventions : Whenever a 
vertex is represented by an empty circle, this means that this vertex neither introduces a weight Wq corresponding to 
its order nor a fugacity g. Consequently those vertices are not counted. Sohd circles correspond to counted vertices 
and therefore introduce factors Wq and g. Combinatorial factors which are due to certain symmetries of the represented 
object will always be displayed exphcitly. Links between vertices will be represented by sohd lineš. As can be seen in 
figuře H the generating function Z will be represented by a bubble, its derivative Z^^^^ by a bubble with a filled circle 
inside and its derivative $ by a bubble with a tail having an empty circle at the end. The tail corresponds to the 
external line of the tree. 

The marked vertex of the ensemble is often called root. Trees generated by Z*^^' are called rooted, those generated 
by $ planted rooted or simply planted. 






O 

(i) (ii) (iii) 

FIG. 3: Graphical representations of the generating functions (a) Z, (b) Z^^^ and (c) $. 

In a similar way one can also define higher derivatives of Z. Each derivative gd/dg introduces a new marked vertex 
and hence another filled circle in the bubble of the graphical representation. Each derivative (1/ g) d / dwi introduces 
a new external line with an empty circle at the end. One could also define derivatives (1/ g) d/dwk which introduce 
an external uncounted vertex connected to the bubble via k links. 

The most fundamental object among all these generating functions is the generating function $ for planted trees. 
With its help one can construct all the others. For example, the combination gwq'^'^ / q\ is a generating function for 
trees with one marked vertex of the order q. If one sums over q one obtains a generating function for trees which have 
just one marked vertex of any order. This is nothing else but Z'^^^ itself. Thus we have : 



^$9^5y(<i>) 



(79) 



q=l 



If one adds a line with an empty end to this marked vertex one obtains the generating function for planted trees. 
The order of the marked vertex to which the line is added consequently increases by one, i.e. q ^ q + 1. Thus the 
corresponding contribution to the sum over q is gWq+i^'^ / q\ : 



(80) 



g=0 



which is a self-consistency equation for $ from which $ can be calculated. Having calculated $ one can insert it to 
( [ŤOl ) and determine Z^-^^ and so on. 






(a) Cb) (c) 

FIG. 4: Graphical representation of gV{<b), gV' {'^) and gV"{'^). 

As we have seen above gV{^) generates trees with one marked vertex, gV'{^) trees with one marked vertex which 
is connected to an external line. One can easily check that the k-th deriváte of V{^) : 



7yW(<I>) = £ 



Wq+k 



$9 



(81) 



9=0 
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generates an ensemble of trees with one marked vertex connected to k external lineš. The graphical representations 
of gV^^^ ($) for fc = 0, 1, 2 is depicted in figuře ^. 

The self-consistency equation (^0|) is illustrated in figuře |^. The content of equation (|7^) emerges automatically 
from the comparison of figuře jŠ] (b) and figuře ^ (a). 





1 
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FIG. 5: Graphical representation of the self-consistency equation (gO|). 

Let US now illustrate how the generating function machinery works solving the classical problém of the tree diagram 
enumeration ||3^. We shall calculate z^r for the case Wt = 1- This is called the Cayley problém. The number of all 
labeled trees with N vertices is given by z^Nl. The self-consistency equation reduces to : 

a> = ge*. (82) 

This equation can be solved for í> : 

N=l 

For the weights Wq = 1 equation ( |79| ) leads to a simple relation between Z^^^ and <i> : 

2« - <i> - .9 - E - E ^-^5^ ■ (84) 

From this relation one can calculate the canonical partition function : 

^7V = ^^ for N>2 (85) 
and the number of labeled tree diagrams to be A^^"^. For large N one can approximate by : 

]\rN-2 



where the last step is due to the comparison with (^9|) . The critical value of the chemical potential and the susceptibility 
exponent take the values — 1 and 7 = 1/2, respectively. It turns out that the value 7 = 1/2 is a generic one. It 
does not change for a wide class of the weights. In the section about the universality classes and singularity types 
two other universality classes of branched polymer models with 7^1/2 will be discussed. 

We will close this appendix with the calculation of the internal two-point function. Similarly as for the partition 
function, it is easier to work with the generating function. Consider tree graphs which are weighted with the fugacity 
g = e~^, and which have two marked vertices separated by n > 1 links. The generating function t/*^^^ defined 
as a sum over all such trees corresponds to the two-point function for the grand-canonical ensemble. Figuře (Iq ) 
illustrates the decomposition of í/^^^ (m? i^to the generating functions gV'{^) and gV"{^) depicted in figuře (H). 
The decomposition is uniquc, since the path connecting the marked vertices is unique. The two bubbles at the ends 
of the chain corrcspond to diagrams of the generating function gV'{^), while the n — 1 ones in between to gV"{^). 
The decomposition leads to the foUowing relation : 

^(2) (^,n) = e-("+i)^ [V ($)]' • [V" ($)]'"~'^ . (87) 

We can also define the internal grand-canonical two-point function for n = 0. In this case the two marked vertices 
lie on top of each other. Thus the two-point function reduces to the one-point function : Q'^^^ il^,n = 0) = Z*^^' = 




Relation ( pl\ ) allows us to find an explicit dependence of the grand-canonical two-point function on n and if we 
íirst solve the self-consistency equation ( ^0|) for We will be rather interested in the scahng behaviour of the 

two-point function near the critical point. 

Let US first illustrate the calculation of the two-point function for the ensemble of trces having the natural weight 
Wt = 1. In this case V'{^) = V"{^) = e*, where is a solution of equation dŠT 



fj. = ^ - log($) . 

In this case the two-point function simplifies to : 

n) = exp {-{n + 1) (a* - ^í>i^l))) = ■ (89) 

The inversion of equation ( ^Šf ) for $ gives <í> = 1 — \/2^/ Afi, when A/i — /i — /íq — /i — 1 0+. Thus the two-point 
function can bc approximatcd in this limit by : 

e(2)(^,n)-exp(-V2(n + l)v/AM) ■ (90) 

We have neglected a term linear in Afj, in the exponent, because we are interested in the limit A/i — > O'^ for which 
A/i ^ \/S/i. What matters in this limit is the leading term in A/i which is related to the large A^-behaviour of the 
underlying canonical-ensemble : 



e(^) (/i, n)^J2 in) e--^ = ^ e''^"^g^^^ {n) e^^^^ , (91) 

N N 



(2) 

where G)^ {n) is the two-point function for the canonical ensemble for trees of size N. In the last formula we split 
into fj, = fj,o + A/í, where /íq is the critical value of /i at which the partition function is singulár. The leading terms 
in A/í of (/íjJT-) are responsible for the scaling behaviour while the next-to-leading ones for finite-size corrections. 

Formula ( |9l] ) is a discrete Laplace transform. Since we are interested in the large iV behaviour of G^^\n), we can 
substitute the discrete by a continuous Laplace transform. The inverse transform then yields : 

g(^^ (n) = e+^"^ -i^ / dA/i íy (^^ gA^AT (92) 

In particular, for the case discussed here (see(|90|)) the exact result rcads : 

g^^^ in) = {2.)-y'e^N-^^'in+l)e^p (-^^) ' (93) 
We inserted /íq = 1 in the last formula. The normalized two-point function can be approximatcd by : 

,(2 



9^\-)- J'.<Í^^ (94) 



The constant c in the last formula is cqual 1. We displayed this constant, because the same formula holds for generic 
trees in generál, but with a constant which depends on the choice of the weights. 

We used two approximations in the last formula. We substituted n -f 1 by n. The difference between the function 
with n and n + í disappears in the large iV- limit. Secondly, the numerator and the denominator in the normalized 

two-point function (M) have a common part (27r)~-'^/'^e^iV~^/^ which does not depend on n. It cancels out. The 

— r . . . (2) 

remaining part : c/N is a normalization constant ensuring g)^ (n) = 1. For finite N there will be some corrections 

to the normalization constant c/N, but these disappear exponentially in the large N limit. 
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Like, for instance, a theoretical physicist performing analytical calculations. 

In the same way one can show that all moments of the two-point function scale uniformly for large A'^ 



X"G);'{X) ~ N 

which means that the two-point function is indeed a function of the scaling argument X/N^^'^ for large A'^. 
[40] For the simulation we used the exact form of G^^' (X) with the constants c and a set equal to one : 



N-1 _kL 

„=i (27rn) 



[41] If all moments of f{x) exist, the size of the centrál region grows faster with n. For example as X* ~ n^^'^ for an even 
distribution. Therefore the gaussian regime sets in faster. 



